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Abstract 

First-principles modeling of a GeO2/Ge(001) interface reveals that sixfold GeC>2, which is derived 
from cristobalite and is different from rutile, dramatically reduces the lattice mismatch at the 
interface and is much more stable than the conventional fourfold interface. Since the grain boundary 
between fourfold and sixfold GeC>2 is unstable, sixfold Ge02 forms a large grain at the interface. 
On the contrary, a comparative study with Si02 demonstrates that Si02 maintains a fourfold 
structure. The sixfold Ge02/Ge interface is shown to be a consequence of the ground-state phase 
of Ge02- In addition, the electronic structure calculation reveals that sixfold Ge02 at the interface 
shifts the valence band maximum far from the interface toward the conduction band. 
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1 With the continued scaling of Si metal-oxide-semiconductor (MOS) devices, it is becoming 

2 increasingly difficult to enhance device performance; therefore, new channel materials have 

3 been explored. Ge is considered as one of the best channel materials for obtaining high 

4 performance MOS devices due to its higher carrier mobility. The narrower band gap of 

5 Ge is also attractive since lower operation voltage lowers energy consumption. To fabricate 
e high performance devices with a Ge channel, one of the most crucial challenges is fabricating 
7 an insulator with superior interface properties since the termination of the surface states is 
s important for device reliability. So far, various insulators have been examined, for example, 
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9 Ge0 2 [l|, l2|, Ge 3 N 4 [3|, and high- k oxides |4|-|6|]. Among these insulators, Ge0 2 is the most 

10 fundamental and important, similar to Si02 in Si MOS technology, because it exists even 

11 in high-A; oxide/Ge interfaces. A considerable number of studies have been conducted on 

12 the process of developing high quality Ge02/Ge interfaces, and several groups have reported 

13 that Ge02/Ge interfaces, fabricated by conventional dry oxidation, have a low interface trap 
u density (mid 10 10 -10 n cm _2 eV _1 ) l|, |2[]. Houssa et al. simulated the density of Ge dangling 

15 bonds at the Ge02/Ge interface as a function of the oxidation temperature, by combining 

16 viscoelastic data of Ge02 and the modified Maxwell's model, and claimed that the density 

17 of Ge dangling bonds is less than that of Si dangling bonds at the Si0 2 /Si interface [7]. 
is Their results are in good agreement with other reported results 

19 Interface stress between a semiconductor and an oxide is considered as one of the origins 

20 of interface defects. Kageshima and Shiraishi predicted that Si atoms are emitted from the 

21 interface to release stress induced by the lattice-constant mismatch between Si02 and Si, 

n 

22 although the dangling bonds remain after Si atom emission [8|. On the other hand, we cal- 

23 culated the probability of Ge atom emission from a Ge02/Ge interface by using the interface 

24 model proposed by Kageshima and Shiraishi and found that hardly any Ge atoms are emit- 

25 ted from the Ge0 2 /Ge interface [9]. We concluded that the flexibility of the O-Ge-0 bonds 

26 contributes to the relaxation of interface stress, resulting in a Ge02/Ge interface that is su- 

27 perior to the Si02/Si one. Watanabe et al. claimed, using the classical molecular-dynamics 

28 simulation, that the narrow equilibrium Ge-O-Ge bond angles contribute to the reduction in 

29 compressive stress in Ge02 films as well as flexible O-Ge-0 bonds lol ]. Although a consider- 

30 able number of first-principles studies on the GeQ 2 /Ge interface have been conducted based 

31 on the calculations on Si02/Si interfaces 0, [ill 15] . the atomic and electronic structures 

32 of the Ge0 2 /Ge interface have not been identified experimentally because Ge0 2 is both 



33 water-soluble and thermally unstable at elevated temperatures. 

34 We propose a complete ordered GeO 2 /Ge(001) interface structure with minimum lattice 

35 mismatch. Our interface model consists of sixfold Ge02, which is derived from cristobalite 

36 and is different from rutile, on a Ge(OOl) substrate. The sixfold Ge02/Ge interface is more 

37 stable by 1.92 eV than the conventional fourfold Ge02/Ge interface. The lattice mismatch 

38 between sixfold Ge02 and Ge (~5%) is much smaller than that between fourfold Ge02 

39 and Ge (~17%). By examining a mixed fourfold and sixfold Ge02/Ge interface model, we 

40 find that sixfold Ge0 2 exists as a large grain at the Ge02/Ge interface. On the contrary, 

41 Si02 at the Si0 2 /Si interface maintains a fourfold structure. A comparative study of the 

42 electronic structures of the sixfold and fourfold Ge02/Ge interfaces shows that the valence 

43 band maximum (VBM) far from the interface varies due to the existence of sixfold Ge02 at 

44 the interface. 

45 Our first-principles calculation method is based on the real-space finite-difference ap- 



46 proach 
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which enables us to determine a self-consistent electronic ground state 



47 with a high degree of accuracy using a timesaving double-grid technique 
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4B norm-conserving pseudopotentials [19] of Troullier and Martins 20[ are used to describe the 



18|. The 



49 electron-ion interaction and are transformed into the computationally efficient Kleinman- 



Bylander separable form 2l|, using the s and p components as nonlocal components for H, 



5i O, and Si, and the s, p, and d components as nonlocal components for Ge. Exchange and 



52 correlation effects are treated using the local density approximation [22] . The coarse grid 

53 spacing of ~0.13 A, which corresponds to the plane wave cutoff energy of ~112 Ry, is used 

54 for all of our calculations. We first examine the atomic structures of Ge0 2 and Si02 bulks in 

55 the cristobalite phases (c-Ge02 and c-Si02) under pressure along the a-axis because these 

56 structures correspond to the directions parallel to the interface when the oxides are piled 

57 on the (001) surface. The c-Ge0 2 (c-Si0 2 ) structure is tetragonal with four Ge0 2 (Si0 2 ) 
5B molecules per unit cell, and 4x4x3 &-point grids in the Brillouin zone are taken into 
59 account. The c-Ge02 structure at the equilibrium point is shown in Fig. Ufa). The equi- 
eo librium lattice parameters of c-Ge02 and c-Si02, which are obtained from first-principles 
ei calculation, are listed in Table HI We then compress c-Ge02 (c-Si02) along the a-axis by 

62 5-25 (5-35)% from the equilibrium lattice constants and optimize the length of the c-axis in 

63 increments of 1% to determine the energy minima. We relax all the atoms until all the force 

64 components drop below 0.05 eV/A. 



65 Figures [5^a) and MJo) show the total energies of c-GeC>2 and c-SiC>2 per molecular unit 

ee as a function of volume. The energy minima of other phases [quartz (□), cristobalite (0)> 

67 and rutile (A)] without any constraints are also depicted for comparison. The zero on the 

68 energy scale is the rutile structure of GeC>2 and the quartz structure of SiC>2. In this figure, 

69 the c-GeO"2 at about 0.78a^ e ° 2 shows a local minimum, where a^ e ° 2 represents the length of 

70 the a-axis of c-Ge02 at the equilibrium point (4.818 A). The atomic structure of the strained 

71 c-Ge02 at the local minimum is shown in Fig. H^b). The c-Ge02 under a certain pressure 

72 transforms into a sixfold structure, which is distinct from the rutile phase, by rotating oxygen 

73 atoms around the Ge atoms. The energy minimum of sixfold Ge02 is lower than that of 

74 fourfold Ge02 since the zero-temperature phase of Ge02 has a sixfold rutile structure. The 

75 critical point between the fourfold and sixfold structures is about 0.85a^°° 2 . The arrow on 

76 the upper horizontal axis corresponds to the lateral length of the Ge(001)-(1 x 1) surface. It 

77 should be noted that Ge02 forms a sixfold structure when the length of the a-axis is equal 
7B to that of the (lxl) surface, while the c-Si02 still maintains a fourfold structure. 

79 We next compare the energetic stability of the sixfold GeO2/Ge(001) interface with the 

so fourfold one since the lattice mismatch between the sixfold Ge02 and Ge(001) surfaces is 

si small (~5%). Figures |3]^a) and Mb) show the fourfold and sixfold Ge02/Ge interfaces, 

82 respectively. The fourfold oxide transforms into a sixfold one by rotating the four oxygen 

83 atoms around one Ge atom in a Ge(001)-(1 x 1) surface unit. The total energy difference 

84 between the fourfold and sixfold Si02/Si interfaces are also calculated for comparison. The 

85 Ge(001)-(v / 2 x \/2) surface is used for the lateral size of the supercell of the Ge02/Ge 
se interface, and the length of the supercell perpendicular to the surface is 5.5a[f e , where a[f e is 
87 the optimized lattice constant of the Ge bulk (5.578 A). The model of the interface includes 
ss seven Ge atomic and two Ge02 molecular layers, and both sides of the surface are simply 

89 terminated with H atoms. Eight /c-points in the lxl lateral unit cell are used for the 

90 Brillouin zone sampling. All the atoms, except the Ge atoms in the bottom-most layer 

91 and the H atoms terminating their dangling bonds, are relaxed. The other computational 

92 details are the same as those used in the bulk calculation. We found that the sixfold Ge02/Ge 

93 interface is more stable by 1.92 eV than the fourfold one because the interface stress between 

94 Ge02 and Ge is released by the phase transition into the dense sixfold structure. On the 

95 other hand, the fourfold c-Si02/Si interface model is preferable by 1.02 eV compared with 

96 the sixfold one. Kageshima and Shiraishi reported that Si atoms at the c-Si0 2 /Si interface 
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97 are emitted to release interface stress, resulting in a quartz-SiC^/Si interface jgj]. On the 

98 other hand, hardly any Ge atoms at the c-GeO^/Ge interface are emitted Our results 

99 indicate that the sixfold structure contributes to the release of interface stress due to the 

100 lattice mismatch between the c-Ge02 and Ge(001) surfaces. 

101 Since the lateral length of the Ge(001)-(1 x 1) surface is longer than that of the sixfold 

102 Ge02 surface but shorter than that of the fourfold c-Ge02 surface in Fig. [2]Ja), there is a 

103 possibility that c-Ge02 on the Ge(001) surface is composed of a mixed fourfold and sixfold 

104 structure. We examine the total energies of the supercell doubling of the Ge(001)-(\/2 x y/2) 

105 surface unit in the two directions, i.e., the supercell contains eight Ge(001)-(1 x 1) units. We 
we respectively replace one and five neighboring (1 x 1) Ge surface units so that 12.5% and 62.5% 
107 of the Ge(001)-(1 x 1) units are composed of the sixfold structures [Figs.|4](a) and IU(b)]. The 
ioB computational procedures are the same as mentioned above. The calculated total energy 
109 differences are summarized in Table [III with respect to the ratio of the sixfold coordination. 

o The fully sixfold Ge02/Ge interface is the most stable, and the mixed interface with the 
12.5% sixfold structure is even more unstable than the fully fourfold Ge02/Ge interface. 
2 The instabilities of the mixed interfaces are attributed to the grain boundaries; the c-axis of 
the fourfold oxidized region is more than 5% longer than that of the sixfold one. This result 
implies that the sixfold oxidized region exists as a large grain at the Ge02/Ge interface. 
Finally, we investigate the effect of the sixfold structure on the variations of the conduction 
e band minimum (CBM) and VBM along the normal direction to the interface. To suppress 

7 the effect of quantum confinement due to the limitation of substrate thickness, a 12-atomic 

8 layer of the Ge(001) substrate is used. The four Ge02 molecular layers are piled on the 
n9 Ge(001) substrate, and other computational details are the same. It is believed that Ge02 

120 is composed of an amorphous structure. Tamura et al. found that the band gap of crystalline 

121 Ge02 is compatible with that of amorphous Ge02 by using the first-principles calculation 



122 23]. Therefore, we calculate the electronic structure with and without a sixfold Ge02 

123 layer inserted between the crystalline fourfold Ge02 layer and Ge(001) substrate. The 

124 grain boundary between the sixfold and fourfold structures parallel to the Ge(001) substrate 

125 is stable, although fivefold Ge atoms exist between the sixfold and fourfold boundaries. 

126 Figures E{a) and ED^b) show the evolution of the CBM and VBM along a direction orthogonal 

127 to the interface plane, respectively. CBM and VBM variations are subtracted from the local 

128 density of states, calculated by integrating them on the plane parallel to the interface based 



129 on p{z,E) = j \if)(r,E)\ 2 dr\\, with a contour of 7.94 xlCT 5 e/eV/A 3 . With the fourfold 

130 Ge02/Ge interface, the VBM is almost complete at about 5 A deep from the interface, 

131 which is similar than that with the SiO^/Si interface [24j, while the CBM saturates within 

132 2 A. Since the band gap of sixfold GeC>2 is narrower than that of fourfold GeC>2 [25], the 

133 valence electrons in the Ge substrate penetrate sixfold Ge02 and the interface dipole emerges. 

134 Therefore, the existence of the sixfold Ge02 layer shifts the VBM far from the interface 

135 toward the conduction band. 



136 In summary, we proposed a sixfold Ge02/Ge interface, in which the lattice mismatch at 

137 the interface is very small (~5%) and which is energetically much more stable than fourfold 
i3B Ge02/Ge interfaces. It should be noted that the sixfold structure was found to be a large 
139 grain at the Ge02/Ge interface after computing the stability of the mixed fourfold and sixfold 
wo Ge02/Ge interface. On the other hand, with Si02, a conventional fourfold structure on the 

141 Si(001) substrate is preferable due to the difficulty in rearranging the rigid O-Si-0 bonds 

142 even in the bulk phase. The electronic structure calculation with and without the sixfold 

143 Ge02 monolayer at the Ge02/Ge interface reveals that the VBM far from the interface in 

144 the ultrathin Ge02 layer (~10 A) depends on the coordination number of Ge02- Our results 

145 provide new insight into a strong candidate of the atomic and electronic structures of the 
we Ge02/Ge interface. We await experimental verification of our prediction. 
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FIG. 1. (Color online) Unit cells of (a) fourfold and (b) sixfold c-Ge02- The solid cube represents 
the unit-cell volume, and the blue (light) and red (dark) circles are Ge and atoms, respectively. 
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FIG. 2. Total energy per molecular unit (m.u.) as a function of volume for (a) c-GeO"2 and (b) 
Si02- Energy minima of other phases [quartz (□), cristobalite (0)> an d rutile (A)] are also shown 
for comparison. The zero on the energy scale is rutile for GeC>2 and quartz for Si02- The upper 
horizontal axes correspond to the lateral lengths of Ge and Si(001)-(1 x 1) surfaces in ao" e ° 2 and 
Og l02 , where a^ e ° 2 and a^ 102 represent the lengths of the a-axes of c-Ge02 and c-Si02 at the 
equilibrium points, respectively. 
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FIG. 3. (Color online) Top views and side views of (a) fourfold and (b) sixfold GeO2/Ge(001) 
interfaces. The blue (light), red (dark), and grey (light small) circles are Ge, O, and H atoms, 
respectively. The dotted square in the top views represents a Ge(001)-(v2 x y/2) surface unit and 
the arrows indicate rotational directions to transform into sixfold structures. 



12 






FIG. 4. (Color online) Top views of mixed fourfold and sixfold structures, (a) 12.5% and (b) 62.5% 
of Ge(001)-(1 x 1) surface units are composed of the sixfold structures. The blue (light) and red 
(dark) circles are Ge and O atoms, respectively. 
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FIG. 5. (Color online) Variations of (a) CBM and (b) VBM along orthogonal direction to interface 
obtained on each atom. Triangles and squares are CBM and VBM of GeO^/Ge interfaces with 
zero and one sixfold Ge02 layers, respectively. Filled (open) symbols represent CBM and VBM on 
Ge (O) atoms. Solid line indicates boundary between Ge02 and Ge(OOl) substrate. 
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TABLES 



TABLE I. Lattice constants of c-Ge02 and c-Si02 at their equilibrium points. The unit is A 



Structure 


a 


c 


c-Ge0 2 


4.818 


7.128 


c-Si0 2 


4.925 


6.828 
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TABLE II. Energy difference between fully fourfold c-GeC>2 structure and various mixing ratios of 
fourfold and sixfold structures with respect to ratio of sixfold structure. All units are in eV. 



Ratio of sixfold structure 


0% 


12.5% 


62.5% 


100% 







0.92 


-0.49 


-7.67 
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